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These notes provide two derivations of the Lorentz-Dirac equation. The first is patterned after 
Landau and Lifshitz and is based on the observation that the half-retarded minus half-advanced 
potential is entirely responsible for the radiation-reaction force. The second is patterned after 
Dirac, and is based upon considerations of energy-momentum conservation; it relies exclusively on 
the retarded potential. The notes conclude with a discussion of the difficulties associated with 
the interpretation of the Lorentz-Dirac equation as an equation of motion for a point charge. The 
presentation is essentially self-contained, but the reader is assumed to possess some elements of 
differential geometry (necessary for the second derivation only). 



I. INTRODUCTION 

The Lorentz-Dirac equation is an equation of motion 
for a charged particle under the influence of an external 
force as well as its own electromagnetic field. The parti- 
cle's world line is described by the relations z a (f), which 
give the particle's coordinates as functions of proper 
time. We let u a (r) = dz a /dr be the four- velocity, and 
a Q (r) = du a /dr is the four-acceleration. The Lorentz- 
Dirac equation is given by 
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(1.1) 



where m is the particle's mass, q its charge, F e " t the 
external force, and a a = da a /dr. Throughout the paper 
we work in relativistic units, so that the speed of light is 
equal to unity. In a Lorentz frame in whi ch the particle 
is momentarily at rest at the time r, Eq. ([^]) reduces to 



q 2 a, 



(1.2) 



where a = da/dt; bold-faced symbols denote the spatial 
components of the corresponding fou r- vectors. The spe- 
cial Lorentz frame in which Eq. (L2) is valid will be re- 
ferred to as the "mom enta rily comoving Lorentz frame" , 
or MCLF. Equation ( |l.2j ) also gives the nonrelativistic 
limit of the Lorentz-Dirac equation. 

My main objective with these notes is to provide a 
self-contained derivation of the Lorentz-Dirac equation. 
In fact, I will present two such derivations, each of which 
relying on a different set of heuristic assumptions. The 
first derivation is patterned after the presentation of Lan- 
dau and Lifshitz Q, which I extend further. The second 
derivation is patterned after Dirac's classic paper 0, al- 
though it differs from it in its technical aspects — my 
own presentation requires much less calculational labour. 
The level of rigour achieved in these notes matches that 
of those two sources. It does not, however, match the 
level achieved in what I consider to be the best reference 
on this topic, the 1980 review article by Teitelboim, Vil- 
larroel, and van Weert ||. Those three references were 
my main sources of inspiration. I have also relied on the 
books by Jackson Q and Rohrlich 



I begin with a review of the fundamental equations 
of electromagnctism in Sec. II. In Sec. Ill I present 
the Landau-Lifshitz derivation of the nonrelativistic 
Lorentz-Dirac equation, Eq. (1.2). This reveals that the 
radiation-reaction part of the vector potential is given by 
|(^4" et — A" dv ), where Af et is the retarded solution to 
the wave equation satisfied by the vector potential, while 
A" dv is the advanced solution. The remaining part of the 
(retarded) potential, ^(^4"^ + ^adv)' does not an?ec t the 
motion of the charged particle. 

Solving the wave equation for Af ct (x) is handled via 
a Green's function whose entire support is on the past 
light cone of the field point x. This introduces a map- 
ping between x and a specific point z on the particle's 
world line, the point at which it intersects x's past light 
cone. The mathematical aspects of this mapping are de- 
veloped in Sec. IV, in preparation for Sec. V, in which 
the potential and field of a point charge are calculated. 
The Landau-Lifshitz calculation of the radiation-reaction 
force is resumed in Sec. VI, this time in a fully relativistic 
context. The calculation is once again based on the half- 
retarded minus half-advanced potential, and it produces 
Eq. (0). _ 

The potential and field at x are naturally expressed 
in terms of u, the particle's proper time when it encoun- 
ters x's past light cone, and r, the distance between x and 
the particle, as measured at that time in the MCLF. This 
motivates the introduction of a (noninertial) coordinate 
system for flat spacetime, based on u, r, and two polar 
angles. This new system is centered on the accelerated 
world line, and Sec. VII is devoted to its explicit con- 
struction. This prepares the way for Dirac's derivation of 
Eq. ({□]), which is presented in Sec. VIII. Dirac's deriva- 
tion is based upon considerations of energy-momentum 
conservation, and it involves only the retarded potential. 

The fact that the Lorentz-Dirac equation involves the 
derivative of the acceleration vector introduces consid- 
erable difficulties in its interpretation as an equation of 
motion. These difficulties are reviewed in Sec. IX; the 
point of view developed there comes largely from the ex- 
cellent discussion contained in Flanagan and Wald [||. 
These difficulties stem from the basic observation that 
point particles cannot be given a fully consistent treat- 



1 



ment in a classical theory of electromagnetism. As a way 
of resolving these difficulties, Landau and Lifshitz take 
the Lorentz-Dirac equation through a reduction-of-order 
procedure, in which q 2 /m is formally treated as a small 
quantity. This procedure is motivated and described in 
detail in Sec. IX. It produces a modified Lorentz-Dirac 
equation, 



(1.3) 



which formally is equivalent to Eq. ( |1 . 1| ) , but is free of 
difficulties. In the MCLF, or in the nonrelativistic limit, 



Eq. (1.3) reduces to 
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(1.4) 



where i^ cx t is the complete (convective) time derivative 
of the external force. 



II. THE FUNDAMENTAL EQUATIONS 

We consider an electromagnetic field F a p produced by 
a point charge q moving in flat spacetime on a world line 
z"(t), where r is proper time. The corresponding current 
density j a is given by 



j a {x) =q / dTu a S(x-z), 



(2.1) 



where u a (r) — dz a /dr is the particle's four- velocity; the 
integration is over the complete world line, and the 5- 
function is four dimensional. Initially we take the world 
line to be arbitrary; our main goal is to find the particle's 
equations of motion. 

The Maxwell field equations are 



(2.2) 



If we express the electromagnetic field in terms of a vector 
potential, 



Fa/3 — Ap^ a — A a fi, 



(2.3) 



and if we adopt the Lorentz gauge, A a a = 0, then the 
field equations take the simple form 



-4irj a , 



(2.4) 



in which □ = i] a/3 d a dp is the wave operator. The 
Minkowski metric diag(— 1, 1, 1, 1) is denoted r) a 0, and 
rj al ^ is its inverse. 

This equation can be solved with the help of a Green's 
function G(x, x'), which satisfies 



UG{x,x') = -4ttS(x - x'); 



the solution is A a (x) = J G(x,x')j°>(x') d 4 x' + A% om (x), 
in which the second term represents a solution to the 



homogeneous equation. The retarded Green's function 
can be expressed as 



G Iet (x,x') = 9(t-t')6{a) 

_ S(t-t' - \x-x'\ 



\x - x'\ 

In the first expression, the quantity a is defined by 
a(x,x') = ^Vapix ~ x') a (x - x') ; 



(2.6) 
(2.7) 

(2.8) 



it is half the spacetime interval between the points x and 
x' . The retarded Green's function has support only on 
the past light cone of the field point x. The second ex- 
pression is obtained from the first by factorizing er; x is 
the spatial projection of the four- vector x a and |x| is its 
magnitude. The advanced Green's function is given by 



G Bdy (x,x')=8tf-t)5(a) 

5(t - t' + \x - x'\) 



(2.9) 
(2.10) 



\x — x \ 

It has support on the future light cone of the field point 



Problem 1. Show that Eqs. (|2.6|)-( |2.10D give the cor- 
rect expressions for the retarded and advanced Green's 
functions. 



The field equations {2.1) imply the charge-conservation 
equation j a a — 0. It is easy to check th at th is equation 
is satisfied by the current density of Eq. (2.1); the calcu- 
lation requires the identity 



u a d a 5(x — z) = — — S(x — z), 

UT 



(2.11) 



which is most easily established in the particle's momen- 
tarily comoving Lorentz frame (MCLF) — the frame in 
which the particle is momentarily at rest at the time r. 
Another relevant conservation equation is 



where 



1.1 V 



(2.12) 



(2.13) 



4tt V M 4 ' 

is the electromagnetic field's stress-energy tensor, while 

= m f dr u a vP S(x - z) (2.14) 



(2.5) 

is the stress-energy tensor of a point particle with mass 



It is easy to show that 
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dpT?i = -F a f = -qj drF a ^6(x - z), (2.15) 
while 



d 3 T; 



a0 



P part 



dra a S(x — z), 



(2.16) 



where a a (r) = du a ' J dr is the p article's ac celer ation. Sub- 
stituting Eq. (fUf]) and (|2.16|) into Eq. (|2.12|) yields the 
Lorentz-force equation, 



(2.17) 



In this equation, F a p must be evaluated at the point 
x = z(t). Sinc e F a p is obviously singular on the world 
line, Eq. (2.17) is not very meaningful as it stands. Our 
aim is therefore to make sense of this equation. 



Problem 2. Show that Eqs. (^2|) and ( |2.17l ) follow 
from the following action principle: 

S = J(~^ F a pF a P + A a jA d 4 x-mj dr. 

Here, j a = q J dX z a S(x-z) and dr = (-g^z 01 ^) 1 / 2 dX, 
where z a = dz a /dX; X is an arbitrary parameter on the 
world line. 



In a specific Lorentz frame, A a can be decomposed 
into a scalar potential $ and a vector potential A. Sim- 
ilarly, F a fj can be decomposed into an electric field E 
and a magnetic field B. In this Lorentz frame, Eq. (2.3) 
becomes 



E = -V$ - 



dA 

~dt' 



B = V x B. 



and Eq. (|2.17| ) reduces to 

ma = q(E + v x B), 



(2.18) 



(2.19) 



where v = dz/dt and a — dv/dt. Also, j a can be de- 
composed into a charge density p and a current density j; 
these ar e ob tained by changing the variable of integration 
in Eq. (|1]) to z°: 

p(t,x) = qS(x - z), j(t,x)=qvS(x-z). (2.20) 

In these expressions, z and v are considered to be func- 
tions of t. 



III. RADIATION REACTION FOR SLOWLY 
MOVING CHARGES 

The assumption that the charge is moving slowly gives 
rise to a clear identification of the radiation-reaction part 
of the electromagnetic field F a p . In this section we will 



show that the part of the vector potential which is re- 
sponsible for the radiation reaction is 



A?Jx) 



1 



(3.1) 



The remaining part, ^(Af ct + A" dv ), must be associated 
with the particle's Coulomb field and does not influence 
the particle's motion. The two parts add up to Af et , 
the correct (retarded) solution to the Maxwell field equa- 
tions. The radiation-reaction potential (|]l]) gives rise to 
a radiation-reaction field F^ , and it is this quantity that 
must be substituted to the right-hand side of Eq. (2.17) 
to obtain the correct equations of motion. Our consider- 
ations in this section will be limited to the nonrelativistic 
limit; the general case will be considered in the following 
sections. 



Using Eqs. (2.7) and (2.10), we find that the retarded 
and advanced solutions to Eq. (2.4) are 



A°{t,x) = 



j a {t-e\x-x'\,x') 



d 3 x', 



(3.2) 



up to the possible addition of a solution to the homoge- 
neous equation; such a term would give rise to an external 
field F^ which can always be introduced at a later stage. 
The parameter e is equal to +1 for the retarded solution, 
and to —1 for the advanced solution. For the time be- 
ing we leave j a arbitrary; it could describe an extended 
charge distribution. We assume that this charge distri- 
bution moves slowly: Let r c be a characteristic length 
scale and t c a characteristic time scale associated with 
the source's motion; we assume v = r c /t c <C 1. We will 
use this condition to simplify the appearance of Eq. fl3.2j). 

Suppose first that we want to evaluate A" far from the 
source. Under the condition r = \x\ ^> \x'\, we have 
| a: — x'\ = r — n ■ x' + 0(r _1 ), where n = x/r is a unit 
vector pointing in the direction of x. If we now expand 
j a in a Taylor series about t — er, treating n ■ x' as a 
small quantity, we obtain 



j a (t e\x x'\,x') = Y/ 1T (n- x') 1 ^ j<*(w, x') 



^ 11 

1=0 



dw l 



where w — t — er represents retarded time if e = 1 and 
advanced time if e = — 1. The slow- motion approxima- 
tion ensures than in this sum, the contribution at order 
I is a factor (r c /t c ) 1 — v l smaller than the leading term 
(I = 0). Substituting this expression for j a into Eq. (3.2) 
yields a multipole expansion for the vector potential: 



i 00 i ii f 
^('>*) = -E7T ni ^ 3 a (w,x>) X ' L d 

1=0 ' J 



(3.3) 



where we use a multi-index L to denote a product of I 
identical factors. For example, ri£ = n ai n a2 ■ ■ -n an and 
summation over a repeated multi-index is understood. 
If e = 1, then rA a is a function of t — r and the vector 
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potential represents an outgoing wave; to such a wave cor- 
responds an outward flux of energy, and outgoing waves 
therefore remove energy from the source. On the other 
hand, if e = — 1 then vA a is a function of t+r and the vec- 
tor potential represents an ingoing wave; to such a wave 
corresponds an inward flux of energy, and ingoing waves 
therefore provide energy to the source. This distinction 
between the retarded (e = 1) and advanced (e = —1) 
solutions is very important for radiation reaction. 



Problem 3. First, show tha t to leading order in a 
slow-motion approximation, Eq. (|3.3|) reduces to 



enp, rA e (t,x)=p, 



where q — J pd 3 x is the total charge and p = J pxd 3 x 
the dipole moment of the charge distribution. The total 
charge is a constant, and p is a function of w = t — 
er; an overdot denotes differentiation with respect to w. 
Second, show that the magnetic field is given by 

rB c = — e n x p, 

and that the electric field satisfies E = e B x n. Third, 
calculate the Poynting vector and show that the rate at 
which energy is flowing out of a sphere of radius r is given 
by 



dE 
dw 



£ 3 P 



Thus, we have an outward flux if e = 1 and an inward 
flux if e = — 1. For a point charge, p = qz and dE j dw = 
e \q 2 a 2 . This is the well-known Larmor formula. 



Suppose now that we want to evaluate A° inside the 
source. In such a situation we may expand j a in a Taylor 
series about t, treating | X X I ELS cL small quantity. This 
gives 



j a (t e\x x'\,x') = £ Ljl \ x x >\l » r{t , x'), 



1=0 



/! 



dt l 



and again we see that the contribution at order I is a 
factor v l smaller than the leading term. In this expression 
we notice that (— e) = +1 if I is even, wh ile (— e) = — e 
if I is odd. Substituting this into Eq. (3.2) gives 



A?(t,x)= ]T ~ />(*,*') I* -ar^dV 

I even 

7 1^1 ** 



I odd 



_> ; T ^-adv)' 



In this expression, the first sum gives |(A' 
and this contribution to A" is the same irrespective of 
the nature of the radiation at infinity, that is, whether 



the waves are outgoing or ingoing. This contribution to 
the vector potential cannot be responsible for the radia- 
tion reaction. On the other hand, the second sum gives 
^(j4" et — A" dv ), and this contribution to A" changes sign 
under a change of boundary conditions at infinity. It is 
clearly this contribution to the vector potential that is 
responsible for the radiation reaction. 

We therefore conclude that t he ra diation-reaction po- 
tential is indeed given by Eq. ( |3 - ip . In a slow-motion 
approximation, this reduces to 



A" T (t, x) 



E 

l odd 



n|r fr(t,x')\x- x '\ i -u 3 x'. 



(3.4) 



The leading contribution to the scalar potential comes 
from I = 3 — the / = 1 term vanishes identically by 
virtue of charge conservation. Thus, 

x) = ~ J p(t, x') \x - x'\ 2 d 3 x'. 

The leading contribution to the vector potential comes 
from 1 = 1: 



A„(t,x) 



d_ 

dt 



j(t,x')d 3 x' 



Both expressions neglect terms which are smaller by fac- 
tors of order v 2 . 

We no w spe cialize to the case of a point charge, and 
use Eqs. fl2.20|). We obtain 



*tr(*,x) - -| |a \X - Z(t) | 2 , A„(t,x) 

Evaluating the time derivatives gives 



d , , 



$rr(t,x) = -{x-z) 



qv ■ a, A„(t, x) 



-qa, 



where a = da/dt. Substituting these expressions into 
Eqs. (2.18), we obtain the radiation-reaction fields, 



E„ 



■ qa, 



B r 



0. 



(3.5) 



Finally, Eq. (|1|) gives us the radiation-reaction force 

2 



F 



9 . 

• q a. 



(3.6) 



The equation of motion of a charged particle is therefore 
ma = F ext + ^q 2 a, in which the first term represents an 
externally applied force. 

We may now substantiate our previous claim that the 
radiation-reaction force is associated with a loss of energy 
in the source. The rate at which the force is doing work 
is given by 



W 



■ v = ^q 2 a ■ v = q 2 



— (a ■ v) — GT 
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Upon averaging, assuming either that the motion is peri- 
odic or unaccelerated at early and late times, we obtain 



(3.7) 



The quantity appearing on the right-hand side is (minus) 
the energy radiated per unit time by the charged particle. 
We therefore have energy balance, on the average. 

Equation (3.6) is the leading term in an expansion of 



the radiation-reaction force in powers of the particle's 
velocity. A c areful examination of the hi gher -order terms 
in Eq. (3.4) reveals that in fact, Eq. ( |3.6j ) is exact in 
the particle's MCLF pj. To generalize to other Lorentz 
frames, we mus t fi nd the unique tensorial equation that 
reduces to Eq. (3.6) in the MCLF. It is easy to check that 
the correct answer is 



(3.8) 



where a a (r) = da a /dr. The equations of motion for a 
point particle of charge q and mass m are therefore 

ma a = F^ t + 2 -q 2 {5 a p + u a u p y i 

where F" xt is the external force. This is the Lorentz- Dirac 
equation. In the following sections we will derive it using 
fully covariant methods. 



Problem 4. Show that Eq. reduces to Eq. fl3.7j ) 

in the momentarily comoving Lorentz frame. Show also 
that the radiation-reaction force can be expressed as 



F° = W{a a - a 2 u a ), where a 2 = a a a° 



IV. LIGHT-CONE MAPPING 

The fully relativistic calculation of the potential and 
field associated with a point charge relies heavily on 
the causal structure of the retarded Green's function 
G Te t(x, x'), which has support on the past light cone of 
the field point x. The field at x therefore depends on 
the state of the charge's motion at the instant at which 
its world line intersects the past light cone. Thus, the 
past light cone defines a natural mapping between the 
field-point x and a specific point z on the world line. In 
this section we develop the mathematical tools associated 
with this mapping. The results derived here will be used 
throughout the following sections. 

Let x be a field point, and z(u) the point at which 
the world line intersects the past light cone of x (see 
Fig. 1). Here, u(x) denotes the specific value of the parti- 
cle's proper time r which corresponds to the intersection 
point; we shall refer to u as the retarded time. Clearly, u 
is determined by a;; it is obtained by solving the equation 




FIG. 1. Light-cone mapping between an arbitrary field 
point x and the point z(u) on the world line. 

= o-(x, u) = \ r, aP [x a - z a (u)} [x? - Z P(u)] , (4.1) 

which states that the points x and z(u) are linked by a 
null geodesic. 

We will need an invariant measure of the distance be- 
tween x and z(u). Consider the scalar quantity 



r(x)^~7 laP [x a -z a (u)] u fi (u). 



(4.2) 



In the MCLF associated with the particle's motion at 
the retarded time, we have r = t — z°(u), which is just 
the time required for light to propagate from z(u) to 
x. Because the speed of light is set to unity, r is also 
the spatial distance between these two points. Thus, the 
invariant r(x) is the distance between x and z(u), as 
measured in the MCLF; we may refer to it as the retarded 
distance between the field point and the particle. Notice 
that because u is determined by x, there is no need to 
make the dependence on u explicit in r: we write r(x) 
and not r{x, u). 

The vector x a — z a (u) is a null vector pointing from 
z(u) to x. It is useful to rescale this vector by a factor 



so as to define a new vector, 



1 



k a (x) = - [x a - z a (u)]. 



By virtue of Eqs. (LI) and (L2), this satisfies 



k a {x)k a {x) = 0, k a {x)u a {u) 



-1; 



(4.3) 



(4.4) 



the second relation provides a convenient normalization 
for the null vector. 

Because x and z{u) are linked by the light-cone map- 
ping, a change of field point x generally comes with a 
change in u. (An exception arises if the displacement is 
directed along the null geodesic linking x to z.) For in- 
finitesimal displacements, it is easy to relate a variation 
in x to the corresponding change in u. Suppose that x 



5 



is displaced to the new field point x + Sx. The new in- 
tersection point in then z(u + Su), and these points are 
related by Eq. (4.1), a(x + Sx, u + Su) = 0. E xpan ding 
t his to first order in Sx and Su, and using Eqs. (4.3) and 
(4.4), we obtain k a Sx a + Su = 0, or 



du 
dx" 



(4.5) 



This result allows us to efficiently differentiate a func- 
tion of x which contains an implicit reference to the 
point u. Let f(x) be such a function. We may make 
the dependence on u explicit by writing f(x) — F(x,u). 
Upon differentiation, we momentarily treat x and u as 
independent variables, and write df = (dF/dx a ) dx a + 
{dF/ du) du. We then use Eq. (4.5) to obtain the full 
derivative with respect to x a : 



df_ 

dx a 



OF 

dx°< 



(4.6) 



This is the differentiation rule under the light-cone map- 
ping. 



As a specific application of Eq. (4.6), it is easy to check 
that the derivative of the retarded distance is given by 



-u a + (1 + rak)k 



(4.7) 



where a& = a a k a is the component of the acceleration 
a a = du a [dr in the direction of k a . It is understood that 
in Eq. (|4.7| ), all world- line quantities (such as u a and a a ) 
are to be evaluated at the retarded time u. Notice that 



Eq. (4.7) implies 



k a (x)r a {x) = 1. 



(4.8) 



As another example, we may use Eq. (L6) to calculate 
the derivative of k a . The result is 



- (g a /3 + kaUp + u a kp - k a kpj - a k 



k a k & . (4.9) 



This implies that k a satisfies the geodesic equation, 
ka^k 13 = 0. We also note that k a a = 2/r. 



Problem 5. Verify Eqs. (fL7|) and (gj). Show that 
r is an affine parameter on the null geodesic linking x to 
z{u). 



discard the possibility of adding a solution to the homo- 
geneous equation, which would give rise to an external 
field F"£.) Substituting Eq. (O), we obtain 



A a {x) = qj d TU a 6{t- z°)S(a), 

where a~(x,z) is defined by Eq. fl2.§D . The (5-function se- 
lects two points on the world line, those which satisfy 
the relation o~{x, z) = 0. One solution gives the retarded 
time u, defined by the condition that z{u) causally pre- 
cedes x; the other gives the advanced time v, defined by 
the condition that x precedes z(v). The (9-function re- 
jects the advanced-time solution, so that A a (x) depends 
on the state of the particle's motion at the retarded time 
only. 

The integral can be evaluated by changing the variable 
of integration to a, writing dr — da /a. We note that 
since a goes from negative to positive values as r passes 
through the retarded time u, a is posi tive a t t — u. In 
fact, a quick calculation using Eqs. (4.1) and (4.2) reveals 
that &(t — u) — r(x). The vector potential is therefore 
given by 



A a (x) 



u a (u) 
' r(x) 



(5.1) 



This is the well-known Lienard- Wichert potential. 

The electromagnetic field F a p is obtained from the vec- 
tor potential by using the differ ent iation rule of Eq. (|4.6j). 
Using u afi = 



k/3 and Eq. ( [4.7|) , we obtain 



2g / 



a k U[ a k{3] 



2q 



«0I> 



(5.2) 



where we have suppressed the dependence of the world- 
line quantities (such as u a , a a , and a k = a a k a ) on the 
retarded time u. The square brackets denote antisym- 
metrization of the indices: a[ a kp\ = ^(a a kp — k a ap). 
(Below, we will use round brackets to denote symmetriza- 
tion of the indices.) It is interesting to note that the part 
of the electromagnetic field which scales as r _1 is propor- 
tional to the particle's acceleration; this is the radiative 
part of the field. On the other hand, the part which scales 
as r~ 2 does not involve the acceleration; this is the bound 
— or Coulomb — part of the field. 



Problem 6. Check the validity of Eq. (|5J). Ver- 
ify that the electromagnetic field satisfies the vacuum 
Maxwell equations, F al3 p = 0, away from the world line. 



V. ELECTROMAGNETIC FIELD OF A 
RELATIVISTIC POINT CHARGE 

The retarded solu tion to Eq. ([Q|), with the current 
density given by Eq. ( 2 . 1 ) , is A a (x) = qj dr u a G rct (x, z), 
in which both z a and u a are functions of proper time r; 
the integration is over the complete world line. (We again 



It is s traightforward to substitute Eq. (|5.2| ) m to 
Eq. (2.13) to calculate the electromagnetic field's stress- 
energy tensor. This computation reveals a natural de- 
composition into radiative and bound components, 



em rad 



a/3 



bnd> 



(5.3) 
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where 



t; 



a0 _ Q 



rad 



a k 



Anr 2 

is the radiative component, and 

,2 



(5.4) 



t; 



a/3 



bnd 



27rr 3 



k^a^ + a k 
2 r 



(k {a i 



0) _ fc" k 13 



2k {a u^ - ^k 13 + - Vaf3 



(5.5) 



is the bound — or Coulomb — component; we use the 
notation a 2 = a a a a . This decomposition is meaning- 
ful because each component is separately conserved away 
from the world line, 



(r ^ 0). (5.6) 



Furthermore, the interpretation of T r "^ as the radiative 
part of the stress-energy tensor is motivated by the fact 
that it scales as r~ 2 and is proportional to k a k^ . 



Problem 7. Check the validity of Eqs. (^-([sTs)) 



Verify that T r "^ is separately conserved away from the 
world line. 



VI. THE RADIATION-REACTION FORCE 

As we have seen in Sec. Ill, the radiation-reaction force 
can be calculated on the basis of a radiation-reaction po- 
tential equal to half the difference between the retarded 
and advanced potentials. The retard ed p otential was cal- 
culated in Sec. V, and given in Eq. (|5.l| ). A similar cal- 
culation reveals that the advanced potential is 



u a {v) 

r a dv(x) ' 



(6.1) 



where v is the advanced time, determined by the relation 
cr(x, z) = and the condition that x precedes z(v ), and 



r a dv{x) = -r] a p[z a (v) - x a ] u p (v) 



(6.2) 



is the advanced distance. 

To construct the radiation-reaction potential, it is ad- 
vantageous to express A2 dv {x) in terms of retarded quan- 
tities. Because we are interested in a field point x very 
close to the world line, this transcription is easily car- 
ried out by Taylor expansion. We shall therefore expand 
At = v — u and r a dv(x) in powers of the retarded dis- 
tance r. It is easy to see that formally, At and r are of 
the same order of smallness. 

We begin by substituting the expansion 



z a (v) = z a + u a At + ~ a a At 2 + ^a a At 3 

+ ±ti a At 4 + 0(At 5 ), 

in which all quantities on the right-hand side are evalu- 
ated at the retarded time u, into the relation u{x, z(v)) = 
0. After some algebra, using the relations a(x, z(u)) = 
and x a — z a (u) = rk a , we obtain 

= 2r At - (1 + ra k ) At 2 - - ra k At 3 

3 
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(rd k +a 2 ) At 4 + (3(At 5 ), 



where, for example, a k = a a k a . This relation implies 
that to leading order, At = 2r — the time delay is twice 
the spatial distance between the world line and the field 
point. Calculating higher-order terms is straightforward, 
if slightly tedious. The result is 



At = 2r 



I-"/./' {a k 2 -^a 2 ~^h k )r 2 + 0(r 3 ) 



(6.3) 

We recall that At stands for v — u, and that all quantities 
on the right-hand side refer to the retarded time u. 

To express r a d v in terms of r, we substitute the pre- 
vious expansion for z a (v) and a similar expansion for 

u a (v), 

u a (v) = u a + a a At + ~ d a At 2 + ^d a At 3 + (9(At 4 ), 



into Eq. (3.2). (Here it is sufficien t to keep terms only 
up to order At 3 .) After using Eq. (6.3), we obtain 



radv = r+-(a 2 + d k )r 3 + 0{r 4 ). 



(6.4) 



We may also use Eq. (3.3) to express u a (v) as an expan- 
sion in powers of r: 

u a {v) = u a + 2a a r + 2(a a - a k a a ) r 2 + 0(r 3 ). (6.5) 

Again, all quantities on the right-hand side refer to the 
retarded time u. 



Problem 8. Check the validity of Eqs. (J6.3|)— (|6.5P . 

Give a physical interpretation to the fact that r ac i v = T 
if the motion is unaccelerated. 



The advance d p otential is ob tained by substituting 
Eqs. (E|) and (E|) into Eq. (Ell). The result is 



A% iv (x)=q— +2qa a + 2q 
r 



a" + a k In 



a — a k a 



r + 0{r 2 ). 
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The radiation-reaction potential is therefore 



A™(x) — — [Af ct (x) — A" dv (x)] 



-qa" 



a a -a k a a -i(a 2 + d^u a 



0(r 2 



(6.6) 



The terms of order r 2 and higher will not contribute to 
the radiation-reaction force, because after differentiation 
they give rise to terms which vanish on the world line. 

Differentiation of the radia tion -reaction potential pro- 
ceeds with the help of Eq. ( |4.6| ) , using such results as 
Eq. (0), ©, and a a .p — ~a a kfj. The calculation is 
straightforward, although it involves a fair amount of al- 
gebra. Many of the terms drop out after antisymmetriza- 
tion of the indices, and we obtain 



x(u,r,B A ) 



FIG. 2. A coordinate system centered on an accelerated 
world line. 



F af3( z ) = -^q^aUf} - U a CLp). 



(6.7) 



Substituting this into the Lorentz-force equation yields 

The equations of motion of a charged particle are there- 
fore 



F 



(6.8) 



where F£^ t is an externally applied force. Thus, the 
Lorentz-Dirac equation follows directly from the realiza- 
tion that the half-retarded minus half-advanced potential 
is responsible for the radiation reaction. The role of the 
remaining part of the retarded potential will be eluci- 
dated in Sec. VIII. 



Problem 9. Check that Eq. (6/7) does i ndee d follow 
from the radiation-reaction potential of Eq. ( |6.6[ ) . Check 
that the Lorentz-Dirac equation is compatible with the 
identity rj a pu a a^ — 0. 



VII. RETARDED COORDINATES 

Our purpose in this section is to develop the mathe- 
matical tools required in an alternative derivation of the 
Lorentz-Dirac equation, to be presented in the follow- 
ing section. We shall construct a coordinate system for 
flat spacetime based on the retarded time u and retarded 
distance r; these coordinates will be centered on the ac- 
celerated world line. A broad class of such coordinate 
systems was considered by Newman and Unti ; ours is 
a specific example. 



The coordinate system (u,r,9 A ), where 8 A = (0,<fr) 
are two polar angles, is constructed as follows. We se- 
lect the point z(u) on the world line, and consider the 
forward light cone of this point (see Fig. 2). To all the 
spacetime events lying on the light cone we assign the 
same coordinate u. The null cone is generated by null 
rays emanating from z{u) and radiating in all possible 
directions. We may choose a specific ray by selecting 
two angles, 6 A , which give the direction of the ray with 
respect to a reference axis. To all the spacetime events 
lying on this null ray we assign the same coordinates 9 A . 
Finally, a specific event on this null ray can be character- 
ized by its affine-parameter distance from the apex; this 
is our fourth coordinate r. 

The point x a (u, r, 6 A ) in spacetime is therefore linked 
to the point z a (u) on the world line by a null geodesic 
of affine-parameter length r; this null ray belongs to the 
light cone labelled by u and is further characterized by 
the angles 6 A specifying its direction on the cone. Let 
the vector k a (u,8 A ) be tangent to this null ray, which 
then admits the description dx a jdr = k a . Because k a is 
a null vector, its normalization is arbitrary, and we are 
free to impose 



Vaff k a {u,9 A )u {u) 



-I; 



(7.1) 



this also determines the normalization of the affine pa- 
rameter r. Its interpretation as the retarded distance 
between x and z(u) comes from the integrated form of 
the geodesic equation, 



x a (u,r,9 A ) = z a (u) + rk a (u,6 A ), 



(7.2) 



which is identical to Eq. (p~3|). Notice, however, the dif- 
ferences in point of view: In Sec. IV, k a was considered 
to be a function of x; here, it is a function of u and 9 A 
which, toget her with r, determine the point x. 

Equation ( [7.2] ) is the desired transformation between 
the inertial system x a and the (nonincrtial) retarded 
coordinates (u, r, 6 A ). So far, the construction is not 
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unique: we have yet to decide how the angles A are to 
be defined. In other words, we must give a prescription 
to locate the polar axis on each of the light cones; the 
angles 9 A will then be the polar angles associated with 
this axis. 

We adopt the following prescription: We consider a 
Lorentz frame (t,x,y,z) which is momentarily comov- 
ing with the particle at the time r = u. In this frame, 
which we shall call the u- frame, u a {u) = S a t , where "=" 
means "equals in the specified Lorentz frame" . We ori- 
ent the spatial axes so that the (purely spatial) accel- 
eration vector points in the direction of the polar axis, 
a a (u) == a(u) S a z , where a(u) is the magnitude of the ac- 
celeration vector. Repeating this construction for every 
point on the world line gives us a polar axis on each of 
the null cones, and the angles 9 A will refer to this choice 
of axis. We let 

k a (u,9 A ) = (l,sin6> cos</>,sin6> sin0,cos6>), (7.3) 

whic h is compatible with the constraints coming from 
Eq. (7.1) and the fact that k a is null. 



Problem 10. Derive the following relations, which 
will be required in Sec. VII: 



— / //" ,/<> : : ll". 

47T 



— / k a k f3 dVL = -g a 

— [ k^k^dCl = - (u a g 
**J 3 V y 



01 j_ uPg^-y + ulg<*P) 



+ 2u a u l3 u~<, 

where dft = sin 9 d9d<fi is the element of solid angle. To 
reduce the work involved, express the null vector as k a — 
u a + n a , where n a is a unit vector orthogonal to u a . Rely 
on symmetry to argue that J n a d£l must vanish, and 
show that J n a n^ d£l must be proportional to g al3 +u a u^. 



calculate the components of k a (u + Su) in the w-frame: 
Using Eq. (|7.3|) , we find k*(u + Su) = 1 + a(u) cos 6* Su, 
k x (u + Su) = s'm9 costp, k y (u + Su) = sin# sin0, and 
k z (u + Su) = cos 9 + a(u) Su. These results imply 

dk a 

— — = (acos9,Q,0,a) = (au u l , 0, 0, a z ). 
ou 

An equivalent tensorial equation is 
dk a 



du 



a k u Q + a a ; 



(7.4) 



it holds in an arbitrary Lorentz frame. With the initial 
data of Eq. (Q, Eq. (JTT^) allows us to find k a (u,9) at 
all times. 



Problem 11. Derive the relations k gdk a /du = and 
u a dk a I du+ak = 0, and show that Eq. ( [7.4[ ) is compatible 
with them. 



We now turn to the calculation of the m etric in the 
coordinates (u,r, A ). Differentiation of Eq. (7^) yields 



dx a =u a du + k a dr + el d9 A , 



where 



and 



dk a 
du 



(1 + ra k )u a + ra a 



dk c 



d9 A ' 



(7.5) 



(7.6) 



(7.7) 



The metric follows by substituting Eq. (7.5) into ds 2 = 
r] a p dx a dx@ and calculating all relevant inner products 
between u a , k a , and e\. For example, the components 
of e\ in the u-frame can easily be obtained from Eq. ( |7.3j ), 
and we infer k a e\ — — u a e\, as well as 



Equation ( |7.3| ) gives an explicit expression for 
k a (u,0 A ) in the u-frame. In the Lorentz frame 
(t' , x', y', z') associated with the time u' , defined so that 
u a (u') = S a t , and a a {u r ) = a(u')S a z ,, k a {u',9 A ) would 
take exactly the same form. This fact can be used to 
obtain a useful differential equation for k a (u, 9 A ), which 
will allow us to work in a single Lorentz frame. To de- 
rive this, we take u' = u + Su and work to first order 
in Su. The relative velocity between the two Lorentz 
frames is u a (u + Su) — u a (u) = a a (u) Su, which allows us 
to deduce the components of u a (u + Su) in the u-frame: 
u a (u + Su) = (1, 0, 0, a Su); we see that the boost param- 
eter is v = a(u) Su. To first order in v, the Lorentz trans- 
formation between the two frames is given by t' ~ t — vz, 
x' = x, y' = y, and z' = z — vt. With this we can 



VapeAe's = diag(r 2 , r 2 sin 2 9); 

these results hold in any Lorentz frame. At the end of a 
straightforward computation, we obtain 



ds 2 = -[{l + ra k ) 2 - 
+ Ira a dud9 / 



2 2 

r a 



du 2 — 2 dudr 
r 2 dn 2 , 



(7- 



where a 2 — a a a a , a k = a a k a , ola — clol^'a, and dfl 2 — 
d9 2 + sin 2 9 dip 2 is the metric on a unit two-sphere. 

When viewed in the retarded coordinates, the met- 
ric of flat spacetime looks horribly complicated — there 
are many nondiagonal terms, and the metric components 
depend on the coordinates. What is then the particular 
advantage of these coordinates? The answer comes from 
Sec. V, in which we saw that the electromagnetic field of 
a charged particle is naturally expressed in terms of u and 







r. This is true as well of the field's stress-energy tensor, 
which will be involved in the calculations of Sec. VIII. 
In fact, these calculations will not involve the metric di- 
rectly. What we shall need, however, is an expression for 
dV , the four-dimensional volume element. In spite of the 
grim appearance of the metric, this is simple: 



dV 



-gd x — r dudrdil, 



(7.9) 



where G?f2 = sin 9 d9d<j) is the element of solid angle. This 
result is established by computing g, the determinant of 
the matrix formed by the metric elements; the square 
root of —g is the Jacobian of the coordinate transforma- 
tion x a (u,t,e A ). 




Problem 12. 

(Q. Show that 
transformation. 



Check the validity of Eqs. ( |7.8| ) and 
i—g is the Jacobian of the coordinate 



Another result we shall need is 



G?£ a = r a r 2 dud£l, 



(7.10) 



which gives the (outward-directed) surface element of 
a three-cylinder r = constant in Minkowski spacetime. 
Here, r a = dr/dx a is the gradient of r in any coordi- 
nate system x a . In the usu al Lorentzian coordinates, r a 
was worked out in Eq. (4.7). In the retarded coordinates, 
r a = S r a . 

The derivation of Eq. ( 7.10 ) proceeds as follows. The 
vectorial surface element on a hypersurface £ is defined 
by dTi a = n a dA, where n a is the hypersurface's unit 
normal (pointing in the outward direction) and dA is the 
three-dimensional surface element. In our application, £ 
is a surface of constant r whose induced metric is ob- 
tained by putting dr = in Eq. (7 



ds| = - [(1 + ra k f - r 2 a 2 ] du 2 



2rctA dud9 J 



dVL 2 



The surface element is then given by dA = y/—gs d 3 x, 
where gs is the determinant of the induced metric. Be- 
cause £ is a surface of constant r, the unit normal must 
be proportional to the gradient of r. We therefore let 
n a = Xr a , and determine A by making sure that n a is 
properly normalized. This gives A~ 2 = g al3 r a rp. In the 
retarded coordinates, A~ 2 = g rr ', a component of the in- 
verse metric. This is given by cofactor(<? rr )/(7, where the 
cofactor of a matrix element is the determinant obtained 
by eliminating the row and column to which this element 
belongs. It is easy to see that cofactor(<7 rr ) = gs, so 
that A = \/g/gj:- Combining these results, we see that 
d£ Q = rig dA reduces to gE q = r ay /g d 3 x, which is just 



Eq. (7.10) 



Problem 13. A particle moving with uniform accel- 
eration in the z direction has a world line described by 



FIG. 3. The world tube E enclosing the world line. The 
dashed lines show a deformed world tube E'. 

z° = a^ 1 sinh(ar), z 1 = z 2 = 0, and z 3 = a -1 cosh(ar), 
where a, the magnitude of the acceleration vector, is a 
constant. These components of z"(r) are given in the 
"laboratory frame" , in which the particle is momentarily 
at rest at t = 0. 

Calculate the components of k a (u,8 A ) in the labora- 
tory frame. Show that the transformation to the retarded 
coordinates is given by 

t = a -1 sinh A + r(cosh A + cos 9 sinh A), 
x = r sin 9 cos <j>, 
y = r sin 9 sin </>, 

z = a -1 cosh A + r(sinh A + cos# cosh A), 

where A = au. Substitute this into the Minkowski metric, 
ds 2 = —dt 2 + dx 2 + dy 2 + dz 2 , and show that you agree 
with Eq. (Q. 



VIII. DIRAC'S DERIVATION 

Dirac's derivation of the radiation-reaction force is 
based upon considerations of energy-momentum conser- 
vation. It goes as follows. We enclose the world line 
within a thin world tube S — a three-cylinder whose 
shape is a priori arbitrary (see Fig. 3) — and we cal- 
culate how much electromagnetic-field momentum flows 
across this surface per unit proper time. We then demand 
that this change in momentum be balanced by a corre- 
sponding change in the particle's momentum, so that the 
total momentum is properly conserved. Provided that the 
particle's momentum is correctly identified, this analysis 
yields the Lorentz-Dirac equation. 

Quite generally, the flux of four-momentum across a 
hypersurface S is given by 



AP a = J T af3 dZf3, 



(8.1) 
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where T a @ is any conserved stress-energy tensor and dT, a 
the outward-directed surface element on E. In our ap- 
plication, E will be the world tube of Fig 3. How does 
AP a depend on the shape of the world tube? The an- 
swer is that it does not: the momentum flowing through 
a deformed tube E' is equal to the flux across the original 
tube E, provided that the tubes begin and terminate on 
the same two-surfaces. (The tubes must therefore have 
the same "caps"; we do not want E' to be longer than E, 
thereby allowing more momentum to go through.) 

The proof of this statement is based on the space- 
time formulation of Gauss' theorem and the fact that 
the stress-energy tensor is conserved. Let AP' a be the 
momentum flux across the deformed tube E'. Let V be 
the four-dimensional region between E' and E; this re- 
gion's boundary is denoted dV and consists of the union 
of E and E'. By Gauss' theorem, we have 



Ap /a _ Apa = I Ta p dE 

'dV 

T ap p dV 

0, 



This is the amount of radiative momentum crossing a 
surface r = constant per unit proper time. In the MCLF, 
this equation reduces to dE/dt = ^q 2 a 2 , the standard 
Larmor formula. 
Going back to Eqs. 



7) and (5.5), we find that the 



radial component of the bound stress-energy tensor is 
given by 



T b„d r /3 



4 7rr 3 



« »i (u a --k a 



Anr 4 



After integration we notice that the most singular terms, 



those which are proportional to r , go away, 
terms, however, give a nonzero result: 



The 



dP a 

du 



2r ' 



(8.3) 



This is the rate of change of the bound momentum. We 
shall have to interpret this result. 



P roblem 14. 

GJ>. 



Check the validity of Eqs. ( |8.2|) and 



which proves the assertion. Because the shape of the 
world tube is irrelevant, we shall make the simplest choice 
and take E to be a hypcrsurface of constant r. The sur- 



face element on this world tube is given by Eq. (7.10) 



We will assume that r is small, so that E lies in the imme- 
diate vicinity of the world line. It will not be necessary, 
however, to take the limit r — > 0. 

We have seen in Sec. V that the stress-energy tensor 
of the electromagnetic field is naturally decomposed into 
radiative and bound parts which are separately c ons erved 
away from the world line. By virtue of Eqs. (4.7) and 
(5.4), the radial component of the radiative stress-energy 
tensor is given by 



The rate of change of elec troma gnet ic momentum is 
given by the sum of Eqs. 



2) and (3.2) 



dP a 
du 



■ q 2 a 2 u° 



(8.4) 



where we have written q 2 /(2r) 



in anticipation of 



our eventual interpretation of the first term. Dirac now 
postulates that the total momentum, P° 



be conserved: 



p mcch. must 



dP 



dP a 

<'iu i mcch q 



du 



du 



(8.5) 



rjiOi/3 „ 

T rad r 



(a 2 -a fc 2 )fc Q , 



47rr 2 

and the flux of radiative momentum is 

,2 



A P a 



4tt 



/ (a 2 — a,k 2 )k a dudQ. 



Notice that the factor r 2 in the surface element has can- 
celed out the factor r~ 2 in the radiative stress-energy 
tensor; the result is independent of r. The rate of mo- 
mentum change is 



dP a 
du 



47T 



J (a 2 -a k 2 )k a dn. 



The angular integration can be handled via the results of 
Problem 10, and we arrive at 



dP; 



rad 



du, 



2 -q 2 a 2 u«. 



(8.2) 



Here, -F^ech * s ^ ne mechanical momentum associated with 



the point particle itself, and Eq. ( p.5[ ) becomes an equa- 
tion of motion if the correct relation between P , and 
the world-line quantities can be identified. This necessi- 
tates a second postulate. 

It would be tempting to make the identification 
J'mech = rriQU a , with mo describing the particle's material 
mass. Combining this with Eqs. (^4|) and ( |3.5[ ), we would 
obtain ma a = —^q 2 a 2 u a , in which m = mo + ra cm would 
be interpreted as the particle's physical mass. This, un- 
fortunately, is a nonsensical result, because the accelera- 
tion cannot be proportional to the four-velocity. 

We must therefore revisit our expression for the me- 
chanical momentum, and allow it to be a more com- 
plicated function of the world-line quantities. We seek 
an expression of the form P™ ech = m u Q + cq 2 a a , in 
which c is a dimensionless constant to be determined. 
(Other possibilities exist, but as Dirac says, "they are 
much more complicated than [this], so that one would 
hardly expect them to apply to a simple thing like an 
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electron".) Combining this with Eqs. (p.4[) and (8.5) 



yields ma a = —q(ca a + |a u a ). Demanding that the 
right-hand side be orthogonal to u a gives c = — § . 

Dirac's second postulate is therefore that the mechan- 
ical momentum must have the form 



-Pmcch = m ou° 



(8.6) 



with mo representing the purely material contribution to 
the particle's mass. This combines with the electromag- 
netic contribution m em = q 2 /(2r) — the electrostatic 
self-energy — to form the particle's physical mass m: 



mo + m em = m + — . 

2r 



(8.7) 



Combining Eqs. (8.4)— (3.7) finally yields the Lorentz- 
Dirac equation, 



■q 2 (a a -a 2 u a ). 



We have seen in Sec. VI that the right-hand side of 
this equation is produced by the half-retarded minus 
half-advanced potential, which is regular on the world 
line. The role of the (singular) half-retarded plus half- 
advanced potential is now clear: it serves to renormalize 
the mass from mo to m; the difference is m em , the elec- 
trostatic self-energy of the particle. 

It should be noted that while this derivation of the 
Lorentz-Dirac equation stays quite close to spirit of 
Dirac's own derivation, it differs from it in its technical 
aspects. Dirac employs a different world-tube construc- 
tion, based on spacelike geodesies (as opposed to null 
geodesies) emanating from the world line. As a result, his 
calculations are much more involved than the ones pre- 
sented here. Another consequence is that Dirac's expres- 
sion for the mechanical momentum dif fers from Eq. ( |3.6[ ) 
it is the expected mou° 



While Eq. (S.6) seems ad hoc 
and strange, it is precisely the expression that results 
from a careful consideration of the "caps" in Fig. 3 - 
it is the mechanical momentum that is appropriate for 
our particular world-tube construction. The proof of this 
statement relies heavily on distribution theory, and is 
presented in Ref. ||. 

The final form of the Lorentz-Dirac equation is ob- 



tained by using the identity a = 



and inserting 



an external- force term on the right-hand side of Eq. ( 3 
This gives 



— p a 

— r ext 







Uf3) 



■-,0 



(8.9) 



If the external force is produced by an external electro- 
magnetic field F^J 3 



, , t , then F? xt = qF"£up. 



IX. DIFFICULTIES OF THE LORENTZ-DIRAC 
EQUATION 

A surprising feature of the Lorentz-Dirac equation is 
that it involves the derivative of the acceleration vector. 
The equations of motion are therefore third-order differ- 
ential equations for z a (r), a very unusual situation that 
necessitates some reflection. For example, an issue that 
must be addressed is the specification of initial data for 
this third-order equation. In the usual case of second- 
order differential equations, the initial data consists of 
the particle's position and velocity at r = 0, and this 
information is sufficient to provide a unique solution. A 
third-order equation requires more, however, and it is not 
clear a priori what the additional piece of initial data 
should be. 

To analyze this problem, let us consider the nonrcl- 
ativistic limit of the Lorentz-Dirac equation, which we 
write in the form 



t a 



1 



where F cxt is an external force and 



2 q 2 



6 m 



(9.1) 



(9.2) 



is a constant with the dimension of time. For the purpose 
of this discussion we take -F G xt to be a given function of 
time. It is easy to check that the general solution to 
Eq. (|Tl) is 



a(t) 



b-^-f F cxt (t')e- t '^dt' 

TOi J-oo 



where we assume that F ext (t) goes to zero in the infinite 
past, sufficiently rapidly that the integral is well defined. 
The constant vector b is not constrained a priori; this 
is the third piece of data that must be specified to com- 
pletely determine the motion of the particle. 

To see how b should be chosen, let us specialize to 
a particularly simple case, in which the external force is 
turned on abruptly at t — and stays constant thereafter: 
F e xt(t) = f6(t), where / is a constant vector. In this 
case we have 



o(i) 



ot/to 



-(' 

m V 



-t/t 



6{t) 



and we see that for arbitrary choices of b, a ~ e*/ to for 
t to- Even though the applied force is constant, the 
acceleration grows exponentially with time. This is the 
problem of runaway solutions, which occurs also in the 
general case. We notice, however, that these unphysical 
solutions can be eliminated if we set b = f/m. Then we 
find 



a(t) = I 
m 



{-t)e 



t/t 



6{t) 
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Although a is now sensibly behaved for t > 0, we see that 
its behaviour is rather strange for t < 0: At a time ~ to 
prior to the time at which the external force switches on, 
the acceleration begins to increase. This is the problem 
of preacceleration, which occurs also in the general case. 



Problem 15. Show that the only non- runaway solu- 
tion to Eq. is 

1 f 00 

a(t) = — / F ext (f + st ) e~ s ds. 
m Jo 

This solution displays preacceleration, because the accel- 
eration at a time t depends on the external force acting 
at later times. 



The problems of runaways and preacceleration cast a 
serious doubt on the validity of the Lorentz-Dirac equa- 
tion. The root of the problem resides with the fact that 
we are trying to describe the motion of a point parti- 
cle within a purely classical theory of electromagnetism. 
This cannot be done consistently. Indeed, a point parti- 
cle cannot be taken too literally in a classical context; it 
must always be considered as an approximation to a non- 
singular, and extended, charge distribution. Essentially, 
the difficulties of the Lorentz-Dirac equation come from 
a neglect to take this observation into account. 

Any extended charge distribution can be characterized 
by its total charge q and a number of higher multipole 
moments. For distances r that are large compared with 
the body's averaged radius £, the electromagnetic field is 
well approximated by the monopole term, q/r 2 , and in 
such circumstances, the point-particle description is ap- 
propriate. For smaller distances, however, this descrip- 
tion loses its usefulness. The point-particle description is 
therefore an approximation in which the internal struc- 
ture of the charge distribution can be considered to be 
irrelevant. This approximation is limited to distances 
r ^> £, and the "world line" of a point charge should 
properly be thought of as an extended world tube whose 
dimensions are nevertheless smaller than any other rele- 
vant scale. Such idealizations are common in physics; in 
hydrodynamics, for example, fluid elements are idealized 
as points. 

We should therefore examine the conditions under 



which Eq. (9.1) is compatible with the restrictions as- 
sociated with a point-particle description. Finite-size 
corrections can be incorporated into the Lorentz-Dirac 
equation, which becomes 



a = — F cxt + t Q a + O(to'&£). 
m 



(9.3) 



This reveals that the original form of the equation will 
be accurate, and the point-particle description valid, if 
ai -C a. To see what this means, let a c ~ f ext /m be a 
characteristic value for the acceleration, and t c a char- 
acteristic time scale over which the acceleration changes. 



Then a ~ a c /t c , a ~ a c /t c 2 , and the correction term will 
be small if 



£< t c 



(9.4) 



Thus, a necessary condition for the validity of the point- 
particle description is that the time scale over which the 
acceleration changes must be much longer than the light- 
travel time across the charge distribution. 

Another condition must be imposed. Any extended 
charge distribution possesses an electrostatic self-energy 
given by m cm ~ q 2 ji. This self-energy contributes to 
the total mass of the charged body: m = mo + m em , 
where mo is the purely material contribution to the mass. 
The classical theory will be consistent if m cm does not 
exceed m: mo must be positive. Our second condition is 
therefore q 2 /£ < m, or 



(9.5) 



Thus, the time constant to must be smaller than the light- 
travel time across t he c harge dist ribution. 

Combining Eqs. (9.4) and ( |9.5| ), we arrive at our most 
important criterion: The point-particle description will 
be valid if and only if 



to < t c 



(9.6) 



Thus, changes in the acceleration must occur over time 
scales that are long compared with to ■ This requirement 
is clearly violated by runaway and preacceleration solu- 
tions. These predi ctio ns therefore lie outside the domain 
of validity of Eq. fl9~T|). 

How do we then go about solving the runaway and 
preacceleration problems? Before answering this, notice 
first that F ext /m ~ a r is the leading term on the right- 
hand side of Eq. ( |9.3|) , that toa is smaller by a factor of 
order (to/t c ), and that the remaining term is smaller still, 
by a factor of order (to/t c ) 2 . The Lorentz-Dirac equation 
can therefore be written as 



a = — F cxt + toa + O{to 2 /t 2 ). 
m 



(9.7) 



Our strategy will be to seek a replacement for this equa- 
tion. (This strategy is the one adopted by Landau and 
Lifshitz; for a thorough review of failed alternatives, 
see Ref. H.) The new equation will be equivalent to 



Eq. (9.7), in the sense that it will differ from it only by 
terms of order (to/t c ) 2 , but it will be free of difficulties. 
In other words, the new equation will be compatible with 
the restrictions associated with a point-particle descrip- 
tion. 

The replacement procedure is sim ple. We start with 
the leading-order version of Eq. (9.7), 



a = —F cxt + O{t /t c ), 
m 



and we differentiate. This gives 



13 



a= — F cxt + 0(t Q /t c 2 ), 
to 



which we substitute back into Eq. (9/7). The result is 



lF cxt + ^F oxt + O(i 2 Ac 2 ), 
m to 



which has the same degree of accuracy as Eq. (p.7\). But 
since this equation no longer involves a, it does not give 
rise to runaway and preacceleration solutions. 

Our conclusion is therefore that the modified Lorentz- 
Dirac equation, 



ma 



toFcxt, 



(9.8) 



is much better suited to govern the motion of a point 
charge, within the restrictions im plied by the point- 
particle description. Equation (9.8) is formally equiva- 
lent to the original Lorentz-Dirac equation, in the sense 
that they are both accurate up to terms of order (to/t c ) 2 . 
But beca use i t is a second-order differential equation for 
z(t), Eq. ( |9.8D does not come with the difficulties associ- 
ated with the original equation. 



T he m ethod by which Eq. (9.7) was transformed into 
Eq. (9.8) is an application of a general technique known 



as reduction of order. This technique can also be applied 
to the relativistic equation, 

a a = — F c ° t + t o (S a + u a u l3 )a' 3 . 

TO f 

To leading order, we have a a = m~ 1 F£ ct , and differen- 
tiation yields d a — m~ 1 F£ x± ^u 13 . Substituting this into 
the original equation, we obtain 



(9.9) 



the relativistic version of the modified Lorentz-Dirac 
equation. 
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Problem 16. 



Show that if the external force is 

a/3 



provided by an external electromagnetic field F c £, then 
the modified Lorentz-Dirac equation takes the form 



ma a = qF"^ t pU 13 + qtc 



+ M5%+u a u )F^ t ^ t y 
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